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ABSTRACT 


There are many important applications of boundary layer flow of the nano fluids over the stretching sheet in manufacturing 
processes. A nanofluid is basically the mixture of nano particles contained in base fluid like water, glycol, and ethylene etc. 
nanofluid helps to improve the fluid’s thermal conductivity and also increases the rate of transfer of heat. In this research 
gyrotactic microorganisms are added along with nano particles within the base fluid. The governing problem flow is tackled 
numerically to explore velocity F’(n), temperature O(n), concentration p(n), motile microorganism density &(n), Nusselt number 
{ Nu} _x, Sherwood number { Sh} _x and local Motile microorganism density parameter { Nn} _x against pertinent physical 
parameters. Figure IV.10 and IV.11 is showing the results with various researches and it is showing the exact agreement with the 
paper results. 


1. INTRODUCTION 

Let us consider a two-dimensional unsteady and 
incompressible flow of fluid due to the stretching sheet. 
Thetransfer of heat of tangent-hyperbolic fluid along the 
wall is with the plane y=0. Fluid flow which is under 
consideration is at y>0. Velocity of stretching sheet is Uw 
and temperature at surface is T,. A magnetic field 
having constant strength By is applied perpendicularly 
to the stretching sheet and magnetic field that is 
induced is so small that it can be neglected. Unsteady 
flow of the fluid and transfer of heat starts at time t=0. 


The geometry of the problem is shown in figure 1.1. 
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Figure 0-1: Geometry of the problem 


2. LITERATURE REVIEW 


There are many important applications of boundary 
layer flow of the nano fluids over the stretching sheet in 
manufacturing processes. A nanofluid is basically the 
mixture of nano particles contained in base fluid like 


water, glycol, and ethylene etc. nanofluid helps to 


— o 
16 International Journal for Modern Trends in Science and Technology 


improve the fluid’s thermal conductivity and also 
increases the rate of transfer of heat. Nanofluids are 
formed because of the dispersion of nano-meter 
substance within the base fluid. Besthapu et. al. [1] 
discussed the effect of viscous dissipation of flow of 
mixed convection which is thermally stratifying MHD 
nanofluid over stretching surface. Sheikholeslami et. al. 
[2] presented the flow of magnetic nano fluid and the 
convective transfer of heat in a porous cavity. 
Sheikholeslami [3] explored the numerical approach for 
transportation of MHD Al203—water nano fluid in 
permeable medium. 

Nowadays bio convective nanofluids over the 
stretching sheet are attracting the researches around the 
world due to its vast industrial applications such as 
biotechnology because they have great mixing property 
and increases the transport of mass. Density gradient of 
motile microorganisms is behind occurrence of 
bioconvection. The microorganisms can swim easily 
and actively in the base fluid. The parameters like 
gravity, oxygen gradient, intensity of light and chemical 
substance are very helpful for microorganisms so that 
they may be able to swim easily in the upward 
Sheikholeslami [4] 


computational approach for entropy and energy 


direction. presented the 
analysis of nanofluid acting under influence of Lorentz 
force acting on a porous medium. Aftab et. al. [5] 
discussed the phase change of Nano confined materials 
for the applications of thermal energy. 

Mosavat et. al.[6] observed study of heat transfer of 
mechanical face seal and fin by using analytical 
method.Lin et. al.[7] introduced the concept of MHD 
unsteady flow of pseudo-plastic nanofluid andtransfer 
of heat in finite thin film with the generation of internal 
heatover stretching surface.Sheikholeslami et. al.[8] 
suggested the novel idea ofmagnetic fieldeffectson the 
unsteady flow of nano fluid and transfer of heat using 
Buongiorno model. Wang et. al. [9] investigatedheat 
transfer and fluid flow of nano particles of copper 
between two parallel porous plates having base fluid of 
Bhattacharyya [10] 
discussed the stagnation-point of boundary layer flow 


water and ethylene-glycol. 
of Casson fluid andtransfer of heat towards a shrinking 
or stretching sheet.Seth et. al.[11] proposed MHD flow 
of stagnation-point andtransfer of heat over 
non-isothermal shrinking or stretching sheet in a porous 


medium source effect.Nadeem et. al.[12] presented a 


numerical solution offlow of non-Newtonian nanofluid 
over stretching sheet.Kumaran et. al.[13] analyzed the 
computational analysis offlows ofMHD Casson and 
Maxwell fluids with cross diffusion. Hamid et. al.[14] 
observed shape effects of nano particles of MoS2 on 
rotating flow of nanofluid along a stretching sheet with 
al.[15] 


suggested analysis ofheat transfer and mass transfer in 


changing thermal conductivity.Hayat et. 
Maxwell fluid in the stagnation region over a stretched 
al.[16] 


Cu-Al203/water hybrid nanofluid in presence of 


sheet.Usmanet. gives the idea of 


nonlinear radiation and thermal 
conductivity.Bhattacharyya et. al.[17] presented transfer 
of reactive solute in MHD boundary layer flow over a 
stretched surface with suction or blowing.Akbar et. 
al.[18] 


inclined asymmetric channel with slip andtransfer of 


analyzed tangent-hyperbolic fluidflow in 
heat.Malik et. al.[19] suggested tangent-hyperbolic 
fluidMHD flow over a stretching cylinder by using the 
method of Keller  box.Akbar et. al.[20] 
proposednumerical solution of boundary layer flow of 
MHDtangent-hyperbolic fluid past a stretching 
sheet.Hayat et. al.[21] presented tangent-hyperbolic 
fluid radiative flow in the presence of convective 
conditions.Usman et. al.[22] discussed the uniform heat 
fluxof Ferro fluids on its flow and heat transferalong a 
smooth plate.Nawaz et. al.[23] proposed study of finite 
element three dimensionalflow of radiative nano 
plasma by using Hall and ion slip currents.Nadeem et. 
al.[24] 
forstagnation flow of non-orthogonal third. order fluid 
al.[25] 


presentedmixed convective flow in inclined square 


gave idea of analysis of nanoparticle 


towards a stretching sheet. Chamkha et. 


filled with nanofluid. Pedley et. al.[26] proposed Bio 
convective flow of a fluid. Kessler et. al.[27] observed 
growth of bioconvective patterns in presence of 
uniformly suspended gyrotactic microorganisms. Hillet. 
al.[28] discussedthe growth of bio convective patterns in 
the presence of gyrotactic microorganisms which are 
suspended in a layer of finite depth. Ghorai et. al.[29] 
proposed the development as well as stability of 
gyrotactic microorganisms in bioconvection. Makinde 
et. al.[30]presented the idea of boundary layer flow of 
nano fluid over a stretching sheet having convective 
al.[31] 
discussed MHDmicropolarnanofluidoverpermeable 

shrinkingsheetin the 


boundary conditions Jayalakshmiet. 


stretching or presence 
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ofNewtonianheating. 


3. METHODOLOGY 


By above consideration, equations of thermal energy, 
total mass, momentum and microorganism can be 


found by nanofluid model as follows: 
V.u = 0, (II. 1) 


+ [p-(1 — BCT -T..))](1 — C) 
— yApN}g, + oJ x B), 


dU ; 
PE = divt + {Cp, 


(I.2) 


= T.L — divq 


+ (22) vr vr] 
T) O Y 


V.j=0, (I. 5) 
In above equations, u = (u, v, 0) is the velocity profile 
of nanofluid, T denotes temperature, C denotes 


concentration, N denotes concentration of 


microorganism, j denotes microorganism flux by 
diffusion, T denotes reference temperature, Dg denotes 
coefficient, Dr denotes 


(V x B) „T.L is 


Brownian diffusion 
thermophoretic diffusion coefficient,J = 
viscous dissipation effect. 
Microorganism flux j is represented by 
j = Nu + NV — D,VN 
Here D, denotes microorganism diffusivity and W, 


denotes as e speed of maximum cell. 
a 
As V= 2, E i 2 and u = (u, v, 0) 
By applying boundary layer assumption, above 


equations becomes 
ðu ðv 
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Boundary conditions 
w=u,, v=0, T=T,, CGPS 
N=N,aty 
=0 (Il. 11)u 
> 0, T > To, C > CE 


N > N.asy > œ (III. 12) 
Using following similarity transformations 


u =-= g= NE zr), Ta Eo 
0n) = Tř- gag a , 
Yi- r (1.13) 


The above equations becomes, 
(1 — n)F" (n) + nWeF"(n)F” 0) — M?F') — F° (n) 
+FM)F'M) -A (ro + ze'o) 
+ [6 — oN, — čRhp]£ 
=0 ae 14) 


e'n) 


TFA —— +10 (n) +Ec(1 -n)F" *(n) += = WeEcPrF" n) 


—A (Zo + xe) — F (men) 


+ FO'N) + NO ME G) + N0 (9) 
=0 (III. 15) 
p'on +e += (rong) Ean) 
=0 (III. 16) 
; (n)nA 
sagon + ER EE) =0 
Boundary conditions becomes 
E(n)=1, Fn)=0, O(n) = 1, 
§(@) =1 as 70 
F(n)>0, 0) > 0, 6) 70, &)>0 as 


yee 


(III. 17) 


om) = 1, 


A. Method of Solution 
The governing equations {III.14} and {III.15} along 
with boundary conditions are highly nonlinear and 


complicated in nature so these are numerically solved 
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by using Runge-kutta Fehlberg method. 
Let 
F =y, F =y, F" = yy 
0 = y3, 0" = yyz, 
Eq. {01.14} and {I1.15}becomes 


(A.1) 
(A. 2) 


n 
(1 = n)yy, +nWeyzyyı + Fyz - y1? — A (yı + 592) 
—M*y,=0 (A.3) 


1 7 f 
592 + Fyz —y18 -4(0 +353) +10 + Ecy, 


n 
+ 3 PrWeEcy,* = (A. 4) 


Boundary conditions are 


F(0) = 0, y,() =1, 
y2(0) = S1, y3 (0) = Sz. 


6(0) = 1. (A. 5) 
(A. 6) 


4. RESULTS AND DISCUSSION 


Graphical results and discussion section is dedicated 


to discover the influence of non-dimensional 
the behavior 


temperature O(n), skin friction Cy and local heat 


parameters on of velocity F'() , 
flux Nuy. Two kinds of physical quantities i-e skin 
friction and reduced Nusselt number are also discussed 
to study the behavior of velocity and temperature 
distributions on the boundaries. 

IV.1lshows the 
distribution F'(ņ) verses magnetic field parameter M to 


Figure behavior of velocity 
distinct values of unsteadiness parameterA. We see in 
the figure that decrease in velocity due to increase inM. 
It is due to the fact that Hartmann number acts as the 
Lorentz force in the field which causes opposition in the 
velocity. So, increase in magnetic field M cause the 
increase in Lorentz force and as a result velocity 
declines. We can see drop in velocity is more forA =0 as 


that for A =1. 


We=0.2. n=0.5 


= tf 4 5 6 7 
Figure 0-1: Velocity distribution F’ (7) verses M to 
different values of A 


Figure IV.2 exhibits the behavior of power-law index 
n on velocity distribution F'(ņ) to distinct values of 
unsteadiness parameterA. In this figure, we see that as 
the values ofn increases, the velocity decreases. This 
behavior is due to Hartmann number M acting as a 
Lorentz force which in turn opposes the velocity 
distribution. Also, for A =0, the decrease in velocity is 


less and for A =1, the decrease in velocity is more. 


We=0.5, M=0.7 


Figure 0-2: Velocity distribution F’(7)versesn to 
distinct values of A 


Figure IV.3 is showing the behavior of velocity 
distribution F'(ņ) verses Weissenberg number We to 
the various values of unsteadiness parameterA. A 
consistent decrease in velocity is seen when the value of 
Weissenberg number We is increased. The physical 
reasoning is due to the drag force which appears 
between the layers of the fluid due to the presence of 
magnetic field. Weissenberg number We is defined as 
the ratio of relaxation time to that of process time. 
Physically, We enhances thickness of the fluid which 
causes in the decline of velocity of the fluid. Also, for A 
=0, velocity decrease is less as that for A =1. 


We = 0, 0.4.0.8 


Figure 0-3: Velocity distribution F’(7)verses We to 
distinct values of A 
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Figure IV.4 bespeaks the consequences of Prandtl 
number Pr to that of temperature distribution (n) to 
distinct values of unsteadiness parameterA. As Prandtl 
number Pr is defined as the ratio of the diffusivity of 
momentum to that of thermal diffusivity. So increasing 
Prandtl number cause the decrease in temperature due 
to the inverse relation of Pr and thermal diffusivity. 
Also for A =0, decrease in temperature is more as that 
for A =1. 


M=0.7. Ec=0.5, n=0.3. We=0.3 


Figure 0-4: Temperature distribution 0(7) versesPr to 
distinct values of A 


Figure IV.5 explores temperature distribution 0(7) 
verses Eckert numberEc to distinct values of heat source 
parameterA. As Ec measures the kinetic energy of fluid 
flow relative to the enthalpy difference. So increasing 
the Ec causes increase in kinetic energy. Physically, heat 
dissipation increases due to increased collision of 
particles of fluid. This dissipated heat is added to the 
system again which causes increase in the temperature. 
Ford = 0, decrease in temperature distribution is more 
and for A = 0.2, decrease in temperature distribution is 


less. 
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Figure 0-5: Temperature distribution 0(n)verses Ec to 
distinct values of A 


Figure IV.6 exhibits the behavior of skin friction 
C,verses Weissenberg number We to the distinct values 
of magnetic field parameter M. Enhancement in values 
of We causes decrease in skin friction because We 
increases the thickness of the fluid which results in the 
decrease in skin friction. Similarly, increasing Hartmann 
number M gives decrease in skin friction. Also, as the 
increases, the value of 


unsteadiness parameter 


coefficient of skin friction increases. 
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Figure 0-6: Coefficient of skin friction Cp verses We to 


distinct values of M 


Figure IV.7 depicts the effect of Weissenberg number 
We on the skin friction Cp to the distinct values of 
power-law indexn. Increasing the Weissenberg number 
We causes the drop in skin friction. For smaller values 
ofn, skin friction is more and for larger values ofn, skin 
friction is less. Also, increasing unsteadiness parameter 


A results in the increasing skin friction. 


kin Friction 
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Figure 0-7: Coefficient of skin friction Cp verses We to 


distinct values of n 


Figure IV.8 shows the influence of Eckert Number Ec 
to the reduced Nusselt number Nu, to the distinct 
values of heat source parameter A. As Eckert number Ec 
is risen-up, reduced Nusselt number is dropped. For the 


increasing values of unsteadiness parameter A gives 
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increase in reduced Nusselt number. 


We=0.3, n=0.2, Pr=6.2, M=0.7 
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Figure 0-8: Nusselt number Nu, verses Ec to distinct 
values of à 
Figure IV.9 bespeaks about the reduced Nusselt 
number Nu, and the power law index n to various 
values of Prandtl numberPr. Reduced Nusselt number 


Nu, increases when Prandtl number Pr is increased. 


Also, enhancing unsteadiness parameter, causes 
enhancement in reduced Nusselt number. 
We=0.3, Ec=0.5, M=0.5 
E1s 5 
É Pr 


n=0, 0.1, 0.2, 0.3 
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Figure 0-9: Nusselt number Nu, verses n to distinct 
values of Pr 


Figure IV.10 is constructed to show the comparison of 
coefficient of skin friction with the existing literature 
when We = 0, n = 0 and for various values ofM. The 
results are showing the excellent agreement with the 


paper results. 
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Figure 0-10: Comparison of present results of 
coefficient of skin friction with previously published 
data when We=0, n=0 and for different values of M 


Figure IV.11 is constructed to show the comparison of 
coefficient of skin friction with the existing literature 
when We=0, 0.3, 0.5 and for the various values of M 
andn. It is observed that increasing the values of M 


gives more accuracy as that of smaller values ofM. 


Review | 


: results | results | i Í results 


Figure 0-11: Comparison of coefficient of skin friction 


with previously published data for different values of 
M, n and We. 


5. CONCLUSION 


A Bio convective flow of a nano fluid past a stretching 
sheet has been investigated in this chapter. Gyrotactic 
microorganisms are added along with nano particles 
within the base fluid. The governing problem flow is 


velocity F'(n) , 
concentration $m) , 


tackled numerically to explore 


temperature O(n) , motile 
microorganism density (n), Nusselt number Nu, , 
Sherwood number Sh, and local motile microorganism 
density parameter Nn, against pertinent physical 
parameters. The main outcomes of this chapter can be 
summarized as: 

e Schmidt 


profileF (q) to dropalong with magnetic parameter M 


number Sc makes the velocity 
but temperature 0(ņ), concentration @(7) and motile 
microorganism density€(4) are enhanced. 

e Velocity profileF (7) decreases by increasing Bio 


convection Lewis parameter Sb while temperature (77) 
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and motile microorganism density (n) are enhanced. 
But there is no visible effect on concentration 
profile (n). 

e Velocity profile F'(ņ) decreases by rising 
Buoyancy ratio parameter Nr but temperature 0 (n) and 
motile microorganism density€() are enhanced. Also, 
there is slight increasing effect on concentration 


profile (n). 


e By raising the values of Thermophoresis 


parameter Nt makes the temperature 0(ņn) , 


concentration (7) and -motile 
densityé(7) to enhance. 


e Increasing the values of Brownian motion 


microorganism 


parameter Nb makes velocity F'(ņ) to decrease but the 
temperature 8(7) and concentration (7) are increased. 
There is no visible effect on motile microorganism 
density€é(7). 

e Increasing bio convection Rayleigh number Rb 
makes temperature@(7), concentration (n) and motile 
microorganism densityé(7) to enhance. 

e Reduced Nusselt number Nu, is increasing for 
Schmidt numberSc and Thermophoresis parameter Nt 
but decreasing for bio convection Rayleigh number Rb, 
Brownian motion parameter Nb, Buoyancy ratio 
parameter Nr and Bio convection Lewis parameter Sb. 

e Sherwood number Sh, and motile 
microorganism density Nn, are showing an increasing 
effect forThermophoresis parameter Nt and decreasing 


effect for other parameters. 


6. FUTURE WORK 


The present research work aims to investigate the flow 
of tangent-hyperbolic fluid past a stretching sheet. 
Unsteady flow and heat transfer was discussed in the 
presence of MHD. Tangent-hyperbolic fluid flow over a 
horizontal stretching sheet was discussed in the 
presence of MHD in the analysis which is the review 
work of Usman et al. An effort was made to extend the 
work done in chapter two by considering the bio 
fluid at the flat 
tangent-hyperbolic fluid. This work may be extended 


convection nano surface of 
by considering the viscous dissipation, velocity slip, 


non-Fourier heat flux and thermal radiation effects 
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